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Abstract. Advanced observational facilities allow to trace back the chemical evolution of the 
Universe, on the one hand, from local objects of different ages and, secondly, by direct ob- 
servations of redshifted objects. The chemical enrichment serves as one of the cornerstones of 
• cosmological evolution. In order to understand this chemical evolution in morphologically differ- 
ent astrophysical objects models are constructed based on analytical descriptions or numerical 
methods. For the comparison of their chemical issues, as there are element abundances, gra- 
^ ■ dients, and ratios, with observations not only the present-day values are used but also their 
temporal evolution from the first era of metal enrichment. Here we will provide some insight 
Pj^ . into basics of chemical evolution models, highlight advancements, and discuss a few applications. 

Keywords, galaxies: abundances, galaxies: evolution, Galaxy: abundances, Galaxy: evolution, 
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o 

Q^) , 1. Introduction 

The evolution of the Universe is clearly manifested by its structure formation and the 
O |. consumption of gas to star formation, the latter leading to the enrichment of chemical 
elements heavier than those stemming from the primordial nuclesynthesis. Knowing four 
facts, the star-formation (SF) rate at any time and its integral over the Hubble time, 
the initial stellar mass function (IMF), the production and yield of each element for each 
particular stellar mass, in each stage of the stellar life, respectively, would well define the 
interstellar enrichment and with some SF delay also the stellar, so that an observational 
t— I ■ trace-back of the chemical evolution (CE) of the Universe would allow to derive one of 
^ ' these parameters definitely and for different specific structures. 

Assuming that a structure like e.g. a galaxy in total or also any galactic region un- 
der consideration behaves like an isolated volume, called a closed box, the longer-living 
£f*^ ■ low-mass stars lead to increase the stellar mass fraction continuously as well as the stel- 
lar remnants' lock-up mass while as a consequence of SF the gas fraction diminishes. 
The basic set of analytical equations for the chemical evolution of galaxies has been 



| formulated by numerous authors (e.g. Talbot & Arnett (1971) Pagel fc Patchett (1975) 
t-h . |Tinsley (1980)||Pagel (1997)|and a lot more since then: see references in the recent papers 



by Prantzos (2008a) and Recchi et al. (2008)| . In this simple model a temporal relation 



for the metallicity as Zi(t) = yi[—ln(ji)] follows where fi — M g (f)/M/ g m is the temporal 
gas mass fraction and yi the yield of element i, i.e. the metallicity release per stellar 
^ . population (see left panel of figH]). 

Unfortunately, nature is not as simple and our present knowledge has a lot of un- 
certainties and even gaps because of models that are reasonably and necessarily much 
simpler than nature and some processes are not yet well understood. A strong simpli- 
fication is implied to the analytical solution by the lifetimes of the stellar factories and 
thus the various delay times of specific element releases. The discrepancy of the so-called 
instantaneous recycling approximation becomes obvious after almost 2-3 Gyrs as demon- 



strated in fig. 12 of Prantzos (2008a) A trace-back of element ratios through long-living 



119 



120 



Gerhard Hensler & Simone Recchi 



simple model 



♦♦♦♦♦♦♦♦^♦♦♦♦♦♦♦♦♦♦♦♦tttliiiliiiit************** 





extreme inflow model 







Figure 1. Evolution in arbitrary time units of an element Zi(t) (green dots), of total mass m 
(blue rhombi), gas mass fractions g (red squares), and stellar mass fraction s (black triangles) 
with respect to t=0 for a closed-box model (left) and for the extreme inflow of primordial gas 
(right) necessary to get the gas mass g constant. 

stars is clearly representing this delay, the most prominent signature being e.g. in the 
[a/Fe]-[Fe/H] trend of the solar vicinity. 

Another simplification can be relaxed by advancing the closed box to open boundaries 
and allowing for gas infall or outflow. Fig. [T] (right panel) demonstrates the extreme case 
of gas infall which is required to balance the gas comsumption. By this, also the stellar 
Z enhancement is reduced by dilution with primordial gas. As a further limitation the 
released gas is assumed to mix instantaneously with the local environment. 

Another effect also enters into the element distribution, namely, the existence of dif- 
ferent gas phases produced by means of various heating processes and the splitting of 
stellar metal injecta to these phases. While the SF is limited to the cool gas alone, cooling 
timescales of dilute hot gas are extremely long, by this, allowing to refrain characteris- 
tic elements from their immediate incorporation into newly formed stars. The ignorance 
concerning such small-scale mixing timescales between the different gas phases until now 
allows only an order-of-magnitude approach for chemical models. Mixing processes are 
not only caused by diffusion in an otherwise static gas but triggered by dynamical effects. 
This means that last but not least basic ingredients to understand and to model the CE 
is gasdynamics. 

Despite of the enormous activities all over the field of CE but because of limited space 
and time in this paper, we cannot cover all aspects and objects under investigation of CE 
and also not refer to the whole bunch of literature, but have to pick out a few papers as 
representative and review the present st ate of modelling the CE of some cosmic structures 
only from the early universe to the local. We therefore apologize for the reasonable 
incompleteness. Starting from the solar vicinity the paper focusses on some hot topics 
as e.g. the galactic halo and its dwarf spheroidal (dSphs) satellite system finishing with 
dwarf gas-rich galaxies. 

Since present-day chemical models are preferentially performed by numerical simula- 
tions a further aspect on the quality and reliability of models should focus on the applied 
numerical methods. Since it is not the aim of this review to also address and discuss the 
validity of numerical methods of galaxy evolution, the interested reader is referred to a 



dedicated recent review by Hensler (2008) 



2. The Milky Way components 

That our Milky Way has experienced such a CE became already obviuos from the find- 
ings by Sandage & Fouts (1987) of the metallicity increase from halo to thin disk stars. 
Although large-scale dynamical effects during the formation and evolution of galaxies, as 
e.g. gas accumulation by infall, galactic winds, etc., are proven by observations, models 
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of the galactic CE are mostly executed for the different galactic components in separa- 
tion only and to some extent implying dynamical processes thru analytical terms, not 
self-consistently but instead tuning them corresponding to a best fit of the observations. 

2.1. The Milky Way disk 

One of the problems within the solar neighbourhood realized at earliest already by Beat- 
rice Tinsley (1974) was the so-called G-dwarf problem, the lack of low-metallicity G 



dwarfs. Although this problem can be simply solved by exponentially declining infall 



(but also outflow or both), Pagel (1987) has already emphasized that additional devia- 



tions from a constant slope in the ln(fi)— Z diagram, the yield y, point towards dynamical 
influences. Moreover, Prantzos (2008a) (his fig. 14) shows that a metal pre-enrichmcnt 
leads to a better fit to the stellar Z — ln(fj,) distribution. 

Therefore, non-dynamical CE models have been conducted with a temporal and radial 



function of a gas infall rate (see e.g. Chiappini et al. (1997) Portinari & Chiosi (1999) ). 



Because also the SF rate is related to the gas surface density S g as ^(r,t) oc vY,*(r,t), 
where v represents a SF efficiency, the infall steers the SF. The exponent k can be set 

relation 1.4. 



to unity or adapted to the original Kennicutt (1989) 



also included, as e.g. the total mass density (see e.g. 



Additional factors are 
Chiappini et al. (1997)[ eq.4). 



The most cited model of this kind by Chiappini et al. (1997) can fit the G-dwarf prob 



lem at best by two infall episodes, a short one with 1 Gyrs timescale only and a long 
one of 8 Gyrs. Astrophysically speaking this so-called two-infall model resembles noth- 
ing exceptional or unexpected than a short initial collapse phase and the general gas 
assembly with long-term decreasing infall rate. Nevertheless, also the SF efficiency as 
another free parameter is changed between both infall eposides. Depending on the stellar 
sample used for the G-dwarf problem the models fit more or less ( |Chiappini et al. 1997[ ), 
but a quantitative comparison of the model abundances at present-day with solar is 
desillusionary and not well surportive for the model since it provides not more than 
a qualitative tendency. The same is discernible for the stellar age-metallicity relation 



(AMR) in Chiappini et al. (1997) 



On short timescales, however, the gas accretion rate should vary locally and affect the 
SF dPflamm-Altenburg fc Hensler 2010] ) and the CE (jHensler et al. 2004aP as well as the 
abundances ( Koeppen & Hensler 200 5| what was also invoked from observations of the 



solar neighbourhood by Knauth et al. (2006) 



How uncertain the deduction of relations can be and to some extent biased by the 
sample taken for studies can also be documented by a comparison of the original AMR of 
stars found by |Twarog (1980j with that derived by Nordstroem et al. (2004)| from a well- 
defined and complete sample of solar vicinity stars. While the traditional results revealed 



the expected tendency but with very strong early enrichment, already Tadross (2003) did 
not find any relation at all from open cluster analyses, while iNordstroem et al. (2004) 



data yield a trend. Showing a large hardly understandable scatter A[Fe/H] by almost 
0.5, even for the average AMR tendency a CE model is hardly constructable without any 
metal pre-enrichment. 

The evolution of radial abundance gradients is reasonable, but values, ranges, and the 
temporal behaviour are also uncertain (see Maciel & Costa, this volume) . Galactic abun- 
dance gradients that are measured from Hn regions as usually by the standard method 
of collisionally excited emission line (CEL, IShaver et al. 1983j) or recently used recom- 
bination lines (RL) in the infrared (JEstcban et al. 2005]) can only yield the present-day 
abundances in the ISM. Over years of determinations the slope has not only been reduced, 
e.g. for log(0/H) from -0.8 per kpc (jShaver et al. 1983|) to -0.44 (jEsteban et al. 2005|) . 
but also a flattening of the gradients in the innermost radii and also outside the solar 
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circle (Vilchcz & Estcban 1996) became available from RLs. From massive stars one can 
derive the same: the chemical abundances of the ISM from which they are born most 
recently. 



While Maciel et al. (2003) could derive a temporal flattening of the oxygen gradient 



over the past 9 Gyrs from -0.11 to -0.06 dex/kpc unsing Planetary Nebulae, from stars 
of different ages Nordstroem et al. (2004) could even attribute a positive slope to disk 
stars older than 10 Gyrs, and again a slight flattening for younger stars from almost -0.1 
dex/kpc to -0.077 since 1.5 Gyrs. 

While the present-day abundance gradients from Hn regions hold for all the determined 
elements, i.e. C, N, O, and a elements, their different radial slopes need explanation. 
Interestingly, the abundance ratio C/O shows also a radial decline what means that C is 
more abundant in O-rich regions. As explanation Carigi et al. (2005) proposed that this 
behaviour is caused by a metal-dependent C yield. On the other hand, also the O yield 
is metal dependent in the sense that at lower metallicities a weaker wind by massive 



stars allow a longer shell-burning production of O (Macder 1992]). Cescutti et al. (2009) 



modelled the chemical evolution of bulge and disk stars with metal-dependent yields by 



Meynet & Maeder (2002b) and could by this also reproduce the enhancement of C/O in 



metal-rich regions like the bulge. 
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Figure 2. Evolution of a 85 Mq star until it s explosion as supernov a typell. Left panel: Wind 
mass loss adapted to evolutionary models by Schaller et al. (1992) Right panel: Time-depen- 
dent abundances of 12 C . 14 N, and 16 Q in the warm Mil region gas after the onset of the WN 
phase at 2.83 Myrs. See|Kroeger et al. (2006) 



For the determination of C and O in Hn regions, another effect has to be studied 
seriously: self-enrichment by Wolf-Rayet stars ( |Kroeger et al. 2 006). At the end of its 
lifetime at t = 3.22 Myr a 85 M star has supplied 0.28 M Q of 14 N, 13.76 M Q of 12 C, 
and 11.12 M Q of 16 0, which are contained in the combined stellar wind bubble/Hn 
region. Since N was released at first in the WN stage it increases slightly in this period 
and is thereafter diluted by the N-poor gas feed. These facts are discernible in the hot gas 
by a first rise and a subsequent decrease of the N abundance ( jKroeger et al. 2006[ ) after 
the transition to the WC stage when C and O are released. The C content in the hot gas 
increases steeply and reaches an overabundance of 38 times solar while the enrichment 
with O is weaker. Due to turbulent mixing of the hot gas with the photo-evaporated shell, 
the warm Hn gas is also enriched but at most for C by not more than a factor of 1.22 
(Fig. [21 right panel). Since the stellar mass range of the Wolf-Rayet phase increases with 
metallicity, the negative slope of the C/O ratio within the galactic disk can be attributed 
to this effect. 
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In a recent study Colavitti et al. (2009) found that the chemical properties of the 



galactic disk are at best reproduced by models with an infall law derived from cos- 
mological but purely Dark Matter simulations. In reality, the disk settling timescale is 
not a simple stellar dynamical one but inherently determined by its thermal and tur- 
bulent energy budget causing a delayed thin-disk formation following the thick disk 
( |Burkert, Truran, fc Hensler 1992] ). 



2.2. Multi-zone models 

As already implied as advancement of simple CE models one can even account for gas ex- 
changes through the boundaries of limited regions under consideration by simply adding 
a corresponding time-dependent (source or sink) term to the gas-mass and metallity 
equations. This was e.g. done in several above-mentioned models with gas infall to the 
galactic disk. Nonetheless, no simple model exists that also accounts for outflow from the 
disk. A first step was the exploration of the effect of the galactic fountain model on radial 



metallicity gradients by Spitoni et al. (2009) The idea is, that the disk gas pushed-up 



by supernovae (SNe) and redistributed while falling back on ballistic paths has to mix 
with the disk gas within the site of impinge and affects the abundance gradient of the 
galactic disk. As general results they found that neither the ballistic delay nor the cooling 
delay of this fountain-driven mixing process alter the metallicity and the radial oxygen 
gradient. 

As well as in this approach also e.g. for infall models, authors preferably avoid the 
complexity of a proper hydrodynamical treatment combined with the chemical and ther- 
mal descriptions by means of an analytical ansatz. From gas infall and expected radial 
gas flows, the above mentioned models mimic the first process by radius-dependent func- 
tions, but drop the latter effect. In order to differentiate between radial disk regions 
and to follow radial flows, multi-zone models were developed e.g. that devide the galactic 
disk into rings and allow for a radial gas flux, are artificially parametrized due to positive 
flows characterizing the inside-out evolution of the disk (jPortinari fc C hiosi 2000]) or are 
consistently determined by the gas pressure gradient. 

Different approaches were performed by coupling structural parts of the Milky Way, e.g. 
halo and disk, and dealing with mass exchange between the regions (jFerrini et al. 1992"! 
1994). By fine-tuning not only the coefficients for mass exchange arbitrarily, but also 
further parameters, as e.g. SF efficiency, etc., rough quantitative agreements with gas 
and stellar mass, abundances, and their radial gradients in the disk could be achieved. 

2.3. Chemo- dynamics 

For all these simple-type CE models it is inherent that their results in comparison with 
observational facts can provide a first insight into amount and timescales on which non- 
local effects have played a role during the evolution. Nonetheless, they are all devoid 
of self-consistency according to the acting processes and the resulting dynamics. Two 
major sources exist, driving the gas dynamics, namely, gravitation and stellar energy 
release. Plasmaphysical processes according to heating and cooling lead to different gas 
phases that have not only different energetics but also behave dynamically differently 
and interact by multiple processes as e.g. by drag forces, shock waves, and turbulence 
as well as energeticly e.g. by heat conduction. Moreover, this means that also SF and 
thermal state of the ISM have to be treated consistently. 

A first-order approach into this direction is provided by multi-phase CE models (also 
implied in Fer rini et al. 19921 11994p where, however, also the interaction terms in a set 
of equations are not self-consistently determined but arbitrarily applied. 

How crucial an appropriate representation of the ISM and its processes is, will be 
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demonstrated in a few of the next sections. Two main ingredients can already be estab- 
lished here: 1) Cool/warm and hot gas phases are coupled energeticly by heat conduction, 
leading to a self- regulation of SF (Kocppcn, Theis, & Hcnsl eT" 1998[ ), local mixture of hot 
and cool gas and, by this, local metal enrichment of cool gas with SN elements. 2) The 
dynamical coupling of both gas phases by drag and mass loading due to evaporation 
hampers the outflow and enhances the cooling of the hot gas. 




s 10 
time [Gyre] 

Figure 3. Metallicity evolution of three different components, cool/warm gas (CM, dark 
hatched), hot intercloud medium (ICM, light hatched), and stars (grey) on average of a Milky 
Way-type chemo- dynamical model (Samland, Hensler & Theis 1997) at 10 kpc galactocentric 
distance within the disk. The right panel is a cut-out of the left one demonstrating how the 
metallicity Z is affected by different processes: 1) Due to infall of primordial gas, Zcm de- 
creases. 2) Zicm increases due to supernova-expelled metals and is always higher than Zcm- 
3) The Zicm decreases therefore, when CM evaporates by heat conduction. 4) Condensation of 
hot gas onto cool clouds increases their Z. Th e stellar Z represe nts an average over stars of all 
ages and is therefore smaller than Zcm- (from Samland (1994) I. 



Up to now as the optimal grid code one can consider the further development of the 



chemo- dynamical scheme by Samland, Hensler & Theis (1997) to 3D and with the stellar 
dynamics for the stars. In addition, a cosmologically growing DM halo is included into the 



simulations by Samland & Gerhard (2003) These models contain all the crucial processes 
of SF self-regulation by stellar feedback, multi-phase ISM, and temporally resolved stellar 
components according to the chemo- dynamical prescription (Hensler 2003). They cannot 
only trace the formation and evolution of the disk galaxies' components but also of 
characteristic chemical abundances and are until now the best self-consistent evolutionary 
models of disk galaxies because the disk formation is included into the global temporal 
galaxy evolution. 

2.4. The Milky Way halo 

Until half a decade ago only poor knowledge on the element abundances of the Milky Way 
halo stars was available. With the wealth of new class 8-10m telescopes, spectroscopy of 
faint and distant halo stars has changed our picture fundamentally. While in the era of 
|Sandage fc Fouts (1987) a halo metallicity [Fe/H] not below -4 could be found and only a 
few below -3, raising doubts on the first population of stars in the halo and their chemical 
fingerprints, recent detections of hyper metal-poor stars (Beers & Christlicb 2005) in the 
galactic halo and their peculiar element abundances (see e.g. IFrebel et al. 2005|) opened 
a wide field of new activities, namely, modelling towards understanding the zero metal- 
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licity nucleosynthesis and studying the formation of the halo. This is necessary in order 
to decide between the two preferred formation scenarios of the Milky Way: the mono- 
lithic collapse model QEggen et al. 1962] ) and the accretion model (jSearle k, Zinn 1978|) . 
A general word of clarification of the two terms should be done here: Even if the accretion 
model has formed an inhomogeneous halo ( |Argast et al. 2000[ |Oey 2000[ ) its mixing of 
infalling gas with stellar ejecta lead to a settling of more or less homogenised gas with 
low angular momentum into the bulge and with larger one into a disk. Due to the higher 
gas column density in the innermost disk due to its smaller area the mostly invoked 
inside-out evolution in the disk is quite plausible. 

Although a dichotomy of halo stars, the inner halo vs. the outer halo fCarollo e t al. 2005)) . 
makes a dissipative formation of the inner halo probable but invokes a dissipationless ac- 
cretion process for the outer halo stars, it is still under debate whether part of the 
present-day dwarf spheroidal galaxies (dSphs) could have served as its building blocks. 
Although most of the dSphs have formed their stellar population in an early short, but 
active SF epoch, the chemical signatures of most of the dSph stars still alive are distinct 
from the stars of each kinematic component of the Milky Way (V enn et al. 2004]) . 

Even the most recent investigation by Prantzos (2008b) which is sold as a clear proof 
that dSphs have served at least partly as building blocks of the Milky Way halo, lacks of 
consistency: SF history, baryonic mass content and metallicity have been assumed from 
the 4 most prominent present-day dSphs, and the accretion timescale was set to half a 
billion years only. 

Since even the CE of the dSphs is not yet conclusively understood, in general, models 
have to be advanced for this morphological galaxy type. 



3. Dwarf galaxies 

3.1. The Milky Way dwarf spheroidal satellite system 

Since the origin and the Dark Matter content of the galactic system of dSph satellites 
is basically questioned by some authors (Mc tz et al. 2009) . general problems are: Can 
low- mass systems survive an intense SF epoch (as expected at the formation)? 

From compilations by Dekel & Woo (2003) and Grebel et al. (2003) dSphs follow a 
mass-metallicity relation. As an extreme case Hensler et al. (2004b) have performed spher- 
ical low-mass galaxies by means of chemo- dynamical simulations in order to study the 
galaxy survival, SF rates, gas loss, and (final) metallicity. They could demonstrate that 
due to the SF self-regulation only short but vehement initial SF epochs occur and lead to 
mass-dependent gas loss. Nonetheless, the dwarf galaxies (DGs) remain gravitationally 
bound with the further issue that more cool gas survives than it is observed, but it forms 
a halo around the visual body. Although the stellar energetic feedback is the driving 
mechanism to expel the gas, its effect is not as dramatic as obtained in semi-analytic 
models (Salva dori et al. 2008P and the amount of unbound mass is considerably lower. 
To get lost, this gas has to be stripped off additionally (Gre bel et al. 2003)) what prob- 
ably happens because of the ram pressure of the galactic halo gas or of tidal stripping 
(jRead et al. 2006]) . Otherwise it can bounce back to the DG and produce subsequent 
events, from a second SF epoch even up to SF oscillations. The external gas reservoir 
around the Scl dSph ([Bouchard et al. 2003|) might whitness this effect. While the models 
make some intermediate-age stars, most stars are made in less than 1 Gyr after forma- 
tion. The dependence of stellar metallicity on the galaxy mass is a quite well confirmed 
observational relation. 

Because of its simplifying one-dimensionality and the neglect of environmental ef- 
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fects like e.g. tidal effects, external gas pressure, gas inflow, etc. more complex models 
are required. The fascinating wealth of data and their precision on stellar ages and 
kinematics, on their chemical abundances and tidal tails of dSphs (for most recent re- 



views see e.g. Koch (2009) and Tolstoy et al. (2009) I have triggered numerous numerical 



models. Although all of them up to date are advanced by 3D hydrodynamics (see e.g. 



Marcolini et al. (2006) and Revaz et al. (2009)1, they still lack of the same aforemen- 



tioned environmental agents. The most comprehensive paper by Revaz et al. (2009) e.g. 
simulates a bunch of DG models with the method of smooth-particle hydronamics (SPH), 
but again all of them in isolation. In their models mainly also sufficient gas mass is re- 
tained and can fuel further SF epochs if it would not be stripped of by ram pressure 
or tidal forces, as the authors mention. Those models that fit the presently best stud- 
ied dSphs Fnx, Car, Scl, and Sex at best, are than assumed as test cases for further 
exploration. Although their results do not deviate too much from the further observa- 
tional data, in addition to the already mentioned neglects, three further caveats must 
be mentioned: 1) If models are selected according to agreement with one or two ob- 
served structural parameters, it is not surprising if also other values would not deviate 
significantly. 2) The numerical mass resolution of the SPH particles is too low to allow 
quantitative issues of galactic winds, heating and cooling, etc. 3) Because of the single 
gas-phase description released metals are too rapidly mixed with the cool gas and the 
metal-enrichment happens too efficiently. Despite these facts, with appropriate initial 
conditions always models in agreement to observations can be found. 

Although the advancement to a two-phase ISM by an SPH code is not trivial and 



inserts various numerical problems, but is not impossible (Berczik, Hensler, et al. 2003 



|Harfst, Theis fe Hensler 2006 Scannapieco et al. 2006 1 , such treatment would be abso- 



lutely necessary in order to achieve reliable results. In addition the chemo- dynamical 
interaction proceeses must be implied. 

3.2. Dwarf Irregular Galaxies 

Dwarf irregular galaxies (dlrrs) are characterized by large gas fractions, often an ongoing 
SF and low metallicities. Until ~ 10 years ago it was believed that at least some of them 
were forming stars now for the very first time. Nowadays it is evident that even the most 
metal-poor ones (like IZwl8) contain stars at least 1 Gyr old ( |Momany et al. 2005[ ). This 
means that SF should have proceeded for at least a few Gyr in dlrrs, albeit at a low 
intensity. The low intensity of the SF in DGs is the best way to explain their chemical 
characteristics, like for instance the low [a/Fe] ratio. As explained in the introduction, the 
[a/Fe] vs. [Fe/H] plot is representative of the different delay in production of a-elcmcnts 
(mostly produced by the short-living massive stars) and the iron (2/3 of which comes 
from longer-living binary systems originating Type la SNe). If the SF duration in a galaxy 
is very short, Type la SNe do not have enough time to restore Fe into the ISM and most 
of the stars will be over-abundant in a-elements compared to Fe. The low average [a/Fe] 
ratios in dlrrs (compared to large galaxies) is a hint of a long-lasting (presumably very 
mild) SF in these galaxies (jlkuta fe Arimoto 2002} ILanfranchi fe Matteucci 2004|) . The 
same trend of [a/Fe] (decreasing as we move towards dwarf galaxies) can be obtained 
by varying the IMF (Recchi et al. 2009]), but the hypothesis of a low SF efficiency in 
dwarf galaxies is important also to explain the downsizing (jCowie et al. 1996)) . namely 
the trend of having older stellar populations as we move towards more massive galaxies, 
indicative of a shorter (but more intense) duration of the SF in large galaxies. 

On the other hand, most people believe that a fundamental role in the chemical evo- 
lution of dlrrs is played by galactic winds. In fact, these systems are characterized by 
shallow potential wells and less energy is required to extract gas from them. Galac- 
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tic winds have the effect of reducing the metallicity of a galaxy in comparison with 
the one predicted by the closed box model. However, detailed numerical simulations 
( |D , Ercole~fc Brighenti 1999; Rcc chi et al. 2006)) have shown that galactic winds are not 
very effective in removing gas from a galaxy. In fact, galactic winds develop vertically 
(along the direction of steepest pressure gradient) whereas the horizontal transport along 
the disk (where most of the gas lie) is very limited. On the other hand, the metals freshly 
produced during a starburst can be easily carried out of the galaxies through a galac- 
tic wind (which will be therefore metal-enhanced) and this, too, contributes to keep the 
metallicity of the galaxy low. Although this overall picture is almost widely accepted, the 
effect of galactic winds depends very sensibly on parameters like galaxy structure and ISM 
properties. In particular, the presence of clouds can hamper the development of galactic 
winds while in the meantime, through evaporation of metal-poor gas, reducing the av- 
erage metallicity of the galaxy. Detailed numerical simulations (Rccchi & Hcnslc r~2007|) 
show that the leakage of metals is not prevented by the presence of clouds and that the 
final metallicity is a few tenths of dex lower than in models without clouds. 

4. Concluding Remarks 

Because of its importance for understanding the evolution of baryonic structures in the 
Universe, their chemical evolution is one of the main focusses of galaxy investigations. 
For simplicity, numerous analytical approaches have been performed for different mor- 
phological units. They can provide a first insight into the temporal and spatial element 
enrichment by stellar nucleosythesis products. 

Since dynamical effects of different gas phases affect the clement enrichment differ- 
entially leading to a redistributions of those elements, but because also plasmaphysical 
processes as e.g. heat conduction and turbulence allow for a mixing of gas, more com- 
plex dynamically and chemically coupled descriptions have to be developed for advanced 
models. The development and application of multi-phase so-called chemo- dynamical de- 
scriptions are on the way. 

How important the implementation of the multi-phase character of the ISM is, can be 
simply understood by the following basic considerations of its inherent nature: In a low- 
mass galaxy the stellar energy release can overcome the binding energy of the whole gas 
reservoir and the gas is totally lost - and by this also the freshly produced metals. Semi- 
analytical (Salvad ori et al. 20 08) or single-phase hydrodynamical (jMacLow fc Ferrara 1999| 
models largely overestimate the gas loss due to the instantaneous energy mixture, but 
to some extent also because of an arbitrarily high SF rate. Although the first-order 
ID chemo- dynamical approach (Hen sler et al. 2004bj) even overestimates the momentum 
transfer by the drag exerted from hot gas flows to the cool clouds and leads to an 
overly more efficient gas removal than in 2D or 3D chemo- dynamical models e.g. of dlrrs 
(jHensler et al. 1999]) . the model description is still closer to reality because it allows the 
outflow of hot gas, while hampered by mass loading and drag of cool infalling intergalactic 
clouds. The amount of expelled metals is reduced (|Recchi fc Hcnsl er~2007j) . 

On massive galaxy scales, as e.g. the evolution of giant ellipticals (gE) ( |Pipino et al. 20 06), 
however, this picture can turn. When hot SNell gas of 10 6 — 10 7 K mixes with a suf- 
ficiently high amount of 10 4 K and cooler gas, a rough estimate gives that with a SF 
efficiency of 10% and a Salpeter IMF about 10% of the stellar mass are returned by 
massive star ejecta, i.e. winds + SNell, what is 1% of the initial SF gas. Accordingly, the 
temperature of the mixed gas then amounts to 2.1 • 10 4 - 1.2 • 10 5 K, i.e. always in the 
range of shortest cooling time, and, by this, remains always bound. A collapsing massive 
galaxy is not stopped and metal-enriched gas is completely following the collapse. Lat- 
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cron, during the evolution this changes due to the gas consumption so that the SN gas 
dominates and gas at above 10 6 K has reasonably to re-expand QPipino et al. 2006]). It is 



not surprising that out of a model grid any model will be found to fit the observations at 
best, but nonetheless the validity of its issued parameters must be questioned. Again the 
simple ID chemo- dynamical model ( | Theis, Burkert, &: Hensler 1992[ ) demonstrates that 
a metallicity of 1-2 times solar can be reached within 1 R e of a gE, while the hot gas 
with 3-4 times solar metallicity extends out to 10 R e . 

The here discussed considerations aim at sensibilizing the reader to the validity of CE 
models. Finally, we wish to mention that the considerations of CE of tidal-tail DGs, of 
the intra-cluster medium, of Ly-alpha structures in the early Universe, and furthermore, 
are also in the focus of crucially important studies and provide further details and to 
some extent initial and boundary conditions for the CE of galaxies. 

Moreover, one should also keep in mind that the stellar yields applied to CE mod- 
els are still uncertain and change due to advanced stellar evolutionary models, tak- 



ing into account stellar rotation (see e.g. Chiappini et al. (2003) Hirschi et al. (2005) 
Meynet et al. (2006)1, binarity (Dc Pon der fc Vanbeveren 2004j) . stellar winds, etc. 
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